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ABSTRACT. In this paper, we derive and analyse a mathematical 
model to assess the impact of three dose vaccination and treament 
as control strategies on the transmission dynamics of rotavirus dis- 
ease. We assume that at least every baby gets the first dose of vac- 
cination and the infected ones are treated. By using the optimal 
control theory, we derive the conditions necessary for optimal con- 
trol of the disease using the Pontryagin’s maximum principle. Nu- 
merical results performed to illustrate our analytical results show 
that in case of an outbreak, with no control infection will disap- 
pear after 75 days, with treatment only, it will take 55 days, with 
vaccination only , it will take between 30 to 40 days and lastly in 
case both treatment and vaccination are implemented the infection 
will disappear within 10 days. Thus both treatment and vaccina- 
tion should be used as controls to fight rotavirus disease since the 
infection takes less time to clear. 


1. INTRODUCTION 


Rotavirus is the most prevalent diarrheal pathogen in young children 
worldwide (Shim et al., 2001). In fact, 95% of the children are infected 
with rotavirus disease by the time they reach age 5 with high incidence 
rate occurring between 4 months and 36 months (Parashar et al., 1998). 
The disease kills an estimated 500,000 children each year and causes 
millions of hospital visits. It is responsible for about 40 percent of all 
diarrheal diseases serious enough to require hospitalizations in young 
children (Ruuska and Vesikari, 1990). The worst burden of rotavirus 
is in Africa, Asia, and Latin America were proper medical care is still 
a problem (CDC, 2009; CDC, 2010 ). 
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Rotavirus spreads through secretions and person to person contact 
(Shim et al., 2001). Contaminated environmental surfaces have also 
been identified as possible mechanism for spreading the disease (Parashar 
et al., 1998). However the disease is curable and preventable. Improve- 
ments in sanitation and hygiene, health education camapigns, treat- 
ment to rehydrate children with severe rotavirus diarrhea can also help 
in the fight of rotavirus disease (Mastretta et al., 2002). 


Given the different control measures available, rotavirus vaccines have 
proved to be the most effective (CDC, 2013). There are two rotavirus 
vaccines on market, that is, RotaTeq and Rotarix (Molholland, 2004). 
Both vaccines have proved to be well much torelated, safe and im- 
munogenic (Vesikari et al., 2004; CDC, 2010). In the United States, 
the introduction of rotavirus vaccines has dramatically reduced seri- 
ous rotavirus infections. In 2008, after the introduction of RotaTeq, 
health care resource utilization for rotavirus disease was reduced by 
almost 90%. Thus, preventing some 50,000 hospital admissions each 
year (CDC, 2010). 


Despite of the available control strategies for the disease, there is an 
urgent need for better understanding of unique parameters in rotavirus 
disease transmission and develop effective optimal strategies for the 
prevention and control of the spread of this disease. 


Alkama et al. (2012) presents an approach that investigates a free ter- 
minal optimal time control of an SIR (Susceptible-Infected-Removed) 
epidemic model. In order to reduce the infected group and increase the 
number of recovered individuals, they present a control simulating vac- 
cination program considering also the minimum duration of a vaccina- 
tion campaign. The optimal control and the optimal final time is found 
using Pontryagin’s maximum principle and the additional transversal- 
ity condition for the terminal time. 


Laarabi et al. (2013) considers a mathematical model of an SIR epi- 
demic model with saturated incidence rate and saturated treatment 
function. They use an optimal vaccination strategies to minimize the 
susceptible and infected individuals as well as maximizing the number 
of recovered individuals. 


Again Tunde and Benyah (2012) considered an SIR model with variable 
size population and formulated an optimal control problem subject to 
the model with vaccination and treatment as controls. Their aim was 
to find the optimal combination of vaccination and treatment strategies 
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that will minimize the cost of the two control measures as well as the 
number of infectives. 


Basing on this studies we formulate an optimal control SIR model with 
vaccination and treatment in the transmission dynamics of rotavirus 
disease. The vaccination in our case is administered three times, that 
is at 2, 4 and 6 months. All doses must be administered not beyond 
8 months. We assume that not all children get all the three booster 
doses however every baby must at least get the first dose. 


1.1. Model Formulation. We develop our model from a SIR basic 
model. The model has seven compartments: M(t) denoting the num- 
ber of breast feeding infants at time t, S(t) denoting the number of 
susceptible children at time t, Vi (t) denoting the number of children 
vaccinated taking one dose between 2 to 8 months, V2(t) denoting the 
population of vaccinated twice between 4 and 8 months, V3(t) denoting 
the population vaccinated thrice, 3rd vaccine between 6 to 8 months, 
I(t) denoting the number of infected children at time t and R(t) denot- 
ing the number of recovered individuals. This is summarized in Table 
1. 


Table 1. Variables of the Model 
































Variables | Description 

M(t) Breast fed infants 

S(t) Susceptible population at time t 

Vi (t) The vaccinated population taking one dose between 2 to 8 months 
Vo(t) The population vaccinated twice between 4 and 8 months 

V3(t) The population vaccinated thrice, 3rd vaccine between 6 to 8 months 
I(t) Infectious population at time t 

R(t) Recovered population at time t 








We consider both breast-fed and non-breast babies. Children are breast- 
fed between 0 — 6 months. The proportion of infants that is breast-fed 
is denoted by q while the birth rate is given by b. The non-breast fed 
fraction enters the susceptible class S(t) at a rate b(1 — q). Maternal 
antibodies wane out at p rate and babies from the breast-feeding class 
enter into susceptible class S(t). 


Let the probability of becoming infected per contact with an infectious 
baby be denoted by 8. The mass infection rate is denoted by STI. 
The daily rate of first dose vaccination between 2 to 8 months is ad- 
ministered at a rate ф. Let т be the daily rate at which children from 
V4 are vaccinated the second time between 4 to 8 months and к is the 
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daily rate at which individuals from V2 are vaccinated for the third 
time between 6 to 8 months. 


Children in the different vaccination groups can come into contact with 
children in infected class, I(t), but the risk between these groups is 
reduced at different rates, 7 as the reduction in risk in case of first dose 
of vaccination 0 < y, € 1 in Vi, т as the reduction in risk for second 
dose of vaccination 0 < 72 < 1 in V3 and lastly, 73 as the reduction in 
risk in case third dose of vaccination 0 < 7з < 1 in V3. Since the more 
times a baby is re-vaccinated the more resistant it becomes, we take 
їз 2122 т. 


Again individuals from /(t) can recover naturally at a rate 0. The 
natural death rate from all the seven classes is denoted by u. In our 
model, we assume that primary vaccination is available but at very low 
levels. Hence we consider both vaccination and treatment campaigns 
as our control strategies. Vaccination campaigns are considered such 
that more children are being vaccinated and the infected ones are being 
treated from the disease. We introduce vaccination as control such that 
the number of susceptible children increase and few become exposed to 
disease. 


We consider u as the fraction of susceptibles being vaccinated per unit 
time for first dose, uz as fraction of being vaccinated for the second dose 
from Vi, из is the fraction of babies vaccinated for the third dose from 
V5. We assume that every child gets the first dose of vaccination. We 
further assume that the vaccine is very effective such that all vaccinated 
children who complete the three doses become immune and they join 
the recovered class R. We again bound our controls with 0 € u; < 0.9, 
for i — 1,2,3 since the entire susceptible population is not vaccinated. 
Children in the infected class are treated at a rate u4 who thereafter 
join the recovery class with 0 € u4 < 1 . In case of no controls, all the 
control parameters are set to zero, that is, ш = из = из = ил = 0. 








'The model flow diagram is described as shown in Figure 1.1 
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FIGURE 1.1. The flow chart describing our model with 
and without controls 


The parameters are summarized in Table 2. 


Table 2. parameters of the Model 

Parameters | Description 

Birth rate 

Death rate 

Proportion of breast feeding infants 

Rate of departure from breast feeding into susceptible class 
Probability of becoming infected per contact 

with an infectious individuals 

The daily rate of first dose vaccination 

The daily rate of second dose of vaccination 

The daily rate of third dose of vaccination 

т Reduction in risk from first vaccine (0 € т € 1) from Vj to I 
n2 Reduction in risk from second vaccine (0 € ro < 1) from V2 to I 
n3 Reduction in risk from third vaccine (0 < 73 < 1) from V3 to I 
0 Natural recovery rate 
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In this model, we derive the best optimal control strategy such that 
the total infection burden is reduced at low minimum cost within a 
period of 100 days thus we have the following set of ordinary differential 
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systems. 
dM 
dS 
p = pM + (1 — q)b — 6S — BSI — uS – ш5 
dV 
т = 08 — MEVI — TV — uV; — uV; 
dV; 
aay ee 
dV: 
— = к\з — 801 — uV; 
dI 
ae BSI + (mVi + Va + n3V3) BL — (0 + u)1 — ual 
dR 
qe rarah | и | us Vi | us V5 


with intial conditions 
M(0) = Me, S(0) = So, VilO) = Vio, (0) = №, V3(0) = 
Үз, I(0)— Io, R(0) = Ro. 


1.2. Optimal Control Analysis. Under this Section, we investigate 
the optimal control efforts needed to control rotavirus disease. Controls 
are represented as functions of time and assigned reasonable upper and 
lower bounds. Vaccination with three booster doses is implemented at 
rate u, (2), u,(t), u4(t) for one dose, two and three doses and u4 is the 
percentage of infected individuals being treated per unit time. Our 
main goal is to minimize the total infectious burden and costs of the 
different controls (vaccination and treatment). The objective func- 
tional J is given as 


OPTIMAL CONTROL OF ROTAVIRUS DISEASE 7 


(1.2) J(u, tg Uz, U4) = 


Т 


(ui, u2, из, ua) Jo 


+ Ba(u4I + u2))dt 


where B,, B,, B,, B, are balancing coefficients or weight constants trans- 
forming the integral into dollars expended over a finite time period of Т 
days. We consider an integrand convex on the set of control parameters 
to indicate non-linear costs potentially arising in case of high treatment 
or vaccination levels. Furthermore we assume unlimited control mea- 
sures and the set of optimal controls (ut (t), už (t), uš (t), u(t)) € Q such 
that 


(L3) — J(uj, uj uj, ui ) = min J (u, (t), u, (t); us (t), u, (t) 
is Lebesgue measurable and is defined as 
(1.4) Q = {(u, (t), u, (t), us (t), ua(t))[0 < u, (t) € t mas 


0 < u(t) <и, 5 0.9,0 < u(t) ш. < 0.9,0 < u(t) <и, 
< Lt € [0,T]} 


< 0.9, 


Since our goal is to characterize an optimal control (и*, už, и*, uš) € Q 
which minimizes the cost of vaccination and the cost of the treatment 
over the specified time interval as well as minimizes the number of 
infectives in a unlimited vaccination scenario. The problem is stated 
as follows: we have to characterize (ut, už, už, иң) € О so that it satisfies 


(1.5) Jul, uz, u, už) = 


T 
in f LE EBS ауа Йй Ge Ci) Bug Анд) 
0 


(ui, u2, ua, u4) 


+ Bs(u4I + uj) }dt 
where О is defined in (1.4) subject to system (1.1). 


Pontryagin’s maximum principle, introduces adjoint functions that al- 
low us to attach our state system i.e. M, S, Vi, V,, Va, I, R differential 
equations (1.1) to our objective functional. The existence of optimal 
controls is guaranteed by standard results in control theory (the inte- 
grand of J is a convex function of the control parameters and the state 
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system satisfies the Lipschitz property with respect to the state vari- 
ables (Fleming and Rishel, 1975; Lenhart and Workman, 2007). The 
necessary conditions that optimal solutions must satisfy are derived 
using Pontryagin’s maximum principle. 


This principle is used to obtain the differential equations for the adjoint 
variables, corresponding boundary conditions and the characterization 
of an optimal control set (ut, и*, и*, už). This characterization gives 
a representation of an optimal control in terms of state and adjoint 
functions. Also this principle converts the problem of minimizing the 
objective functional subject to the state system into minimizing the 
Hamiltonian with respect to the controls (bounded measurable func- 


tions) at each time t. 


According to (1.1) and (1.2) the Hamiltonian for our optimal control 
problem is: 


(1.6) 
H =I + Bo(u48 + u?) + В, (и + из) + Ba(u3 Va + uz) + Ba(u4I + u$)- 
Au(t)(ba — pM (t) — uM (t))+ 
As(t)(pM + (1 — q)b — 6S — BSI — uS —uiS)— 


Av, ()(68 — mBV,I ту — pV, — uV) + 

Av, (D (ти — n28V,I — u,V, — KV, — nV;)-- 

Av, (0) (6% — 13/8131 — рУз) + 

Ar(t)(BST + (mVi + 1 + m3V3)81 — (0 + и)1 — и41)+ 








An(t)(01 + u4I — uR | и | ТРАГА | ua V3). 


where Ам, Às, Av,;Av,,Av,,Ar, Ав are the adjoint functions associated 
with their respective states. Note that, in H, each adjoint function 
multiplies the right-hand side of the differential equation of its corre- 
sponding state function. The first term in H comes from the integrand 
of the objective functional. 


Given an optimal control set (už, už, uš, u?) and corresponding states 
(M*,S*, V*, VŽ, VŽ, I*, R*), there exist adjoint functions 
AM; Às, Ау, › rv, ; XV,» Аг, Ав satisfying 


dàm _ ӨН dis ӘН d ӧн dM, ӘН dày, _ 











dt OM’ d OS’ dt OW dt OV,’ dt 
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4м _ ӘН dg — OH dx, ӘН dix, ӘН Ax, _ 
dt Of’ dt OR dt OX,’ dt — OXy dt 
_ OH 
OX3 
That is 
dA, _ 
dt =. Ам (p T Hu) Asp 
dX 
A ER Às (Ф + 81(t) tA- ui (t)) ЗЕ А, Ф = A,BI(t) un Ви, (t) А Акма (2) 
ал 
m = rv, (mG (t) тт ит Uy, (t)) > Ay, T E Arm BI(t) m Bu, (t) m Agu2 
ал 
T ж А, (n,BI(t) RFH u,(t)) = AK = Xn, BI(t) Ts B,u,(t) a= Лиз 
dà, 
л = Ау, (7,61 (t) Y и) e A DI(t) 
dX 
Tt = ASBS(t) + dy, n, BV, (t) + Ay, MB V(t) + Av, aB V: (0) ee 
ACCS CE) + m, V, (t) + m, V, (t) + m, V;(t)8 — 0 — us — u) — A&(0 + Ua) 
— Взид 
dX 
apt 


The transversality conditions i.e. the final time boundary conditions for 
the adjoint variables Ам, As, Av, , Ау,, Àv, , Ат, An are zero since there is 
no dependence on the states at the final time in the objective functional. 


The Hamiltonian is minimized with respect to the control (at the 
optimal control set) thus we differentiate H with respect to u,, u,, Us, и, 
respectively in the interior of О to obtain the optimality conditions that 


10 NAMAWEJJE, GHOSH, FERRARI, LUBOOBI 




















follows: 
OH 
Du, = В,(5 + 2u,) —A,S+A,5 = 0 
Н 
- = В, (V + 2u,) — A, V, +А, И, = 0 
ди, 1 
(1.8) aH 
ðu, = B, (V, T 2u, ) ui ^v, Vo v Ar V2 =0 
oH 
ди; = B (I + 2u,)— AM FAI = 0 


A 


Solving for u, = U,,u, = û,, U, = à,,u, = U, from the above set of 
equations, the optimality conditions are obtained as follows: 














бый 
: 2B, 
Ama bd (Av, — àg — BM 
zii 2B, 
1.9 
( пе (Ар = Àn =. By)V2 
B 2Bo 
ОАЕ 
5 2B3 


Using the standard argument for control bounds, we arrive at the fol- 
lowing expression for the optimal control function: 


u*(t) = min{maxr{0, i, (©) }, и, } 
u(t) = min([maz10,à,(t)),u,...) 
(1.10) u;(t) = тїп{тахт{0,@„(Ё)}, us} 
už (t) = min(maz(0, à,(t)),u,....) 


From (1.10) the following optimality system is obtained. 
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dM 
= М М 

di bq— p uM, 
ds | À 
Ge ^M + (— q)b — 68 — BST — u$ — min{mar{0, à, (£)],u, pas} 

dV, : Я 
us oS — s 8ViI — TV, — pV — min{mazx{0, û, (t) }, и, И, 
ау Р : 
2 = тїї mÉV,I — min(max(0, à, (0), ta „a, V, — &V, — nV, 
dV: 
Gp = к — ns BVal — Vs 
E І х 
— = BSI + (mV: + т + nsVs)BI — (0 + p)! — тіт{тал:{0, à, (t)}, ч, , 
dR f А ; 2 
es ӨТ + min(maz(0,à,(t)],u,,— 1 — HR + min([maz(0,à,(t)),u,.— 1S 
min{max{0, à,(t)),u, Vi + min(maz(0, à,(t)),u, .— Va 
(1.11) 
dX, 
"di. pe Ам (p RE Hu) Asp 
dA; К 

dt = 5(Ф + BI(t) tnu- min{mazx{0, и, (t)}, а. 0..7) EX ^v ф m A, BI(t)— 
B,min(maz(0, à, (t) ),u,.— } = Aymin(maz(0,à,(t)),u,,. } 
dày, 

"x =e ‚ (m 8I(t) +7 + p+ miimaz(0,ü,(0)),u, }) — Ау, TO 
А810) — B, min{max{0,,(t)}, us} —rqmin{mar{0,@,()}.u,,_} 
dày, | А 

T = ày, (BI (t) + & + u + min{max{O, û, (t)}, us, H Av, k— 

Am 81(t) — B, min(maz(0, à,(0)). u,...  — Атїп{таа{0,й,(@}, us} 
dà; 

ps — ^v. (n, BIC) Tm и) = A ns bI (t) 

ал, 

t = A; S(t) + ^v n BV, (t) F Av, т ВУ, (t) + rv, ПЗУ» (t) m m 
A, (S(t) + m, V, (t) + n, V; (0) + n V, (t))8 — 0 — min{maz{0, а, (t)}, ts nas} — H) 
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Ак(0 + min(maz(0,à,(t)),u,.— }) — Bsmin(maz(0,6,(t)),u,.— } 
ал, 
dp TORNE 


with initial conditions 


M(0) = Mo, 5(0) = So, №(0) = Vio, (0) = Иж, V3(0) = 
V30, I(0) = Io, R(0) = Ко 


and А„(Т) = 0, ^„(Т) =0, à, (T) = 0, Ay, (T) = 0, Aå. (T) = 
0, А„(Т)=0, ^„(Т)=0. 

Now, the state system of differential equations and the adjoint sys- 
tem of differential equations together with the control characterization 
above form the optimality system to be solved numerically. Since the 
state equations have initial conditions and the adjoint equations have 
final time conditions, we cannot solve the optimality system directly by 
only sweeping forward in time. Thus an iterative algorithm “ forward- 
backward sweep method ” is used. An initial estimate for the control 
is made. The state system is then solved forward in time from the 
dynamics using RK4 method. Resulting state values are placed in the 
right hand side of the adjoint differential equations. Then the adjoint 
system with given final conditions is solved backward in time, again 
employing RK4. Both state and adjoint values are used to update the 
control using characterization and then the process is repeated. This 
iterative process terminates when current state, adjoint and control 
values converge sufficiently. 


2. NUMERICAL ANALYSIS OF THE MODEL 


In this Section, we study the optimal control model (1.1) numeri- 
cally. We solve the model using the forward backward sweep method 
(Lenhart and Workman, 2007). Using the iterative method, we solve 
the optimality system derived in (1.11). We examine the model with 
three doses of vaccination controls, ш, шо, из and treatment u4 on 
the spread of rotavirus. In our numerical analysis, we investigate and 
compare numerical results with different scenarios, and these include: 


(i) When treatment is optimized and vaccination is set to zero, 
that is, u; = ио = из = 0 

(ii) When controls u1, из, из, that is, vaccination is optimized while 
treatment is set to zero. 
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(iii) When both controls are optimized, that is, when both treatment 
and vaccination are optimized. 
(iv) Effect of no control, with treatment only, with vaccination only, 
and both controls on the infected class. 
We further assume the weight factors for vaccination is the same for 
all doses and treatment is slightly lower. The cost factors associated 
with treatment include, oral rehydration drugs, medical examinations 


and hospitalization. The parameter values used in the simulation of 
the model are stated in Table 3a and Table 3b. 


Table 3a. State variables and parameters of the Model 









































State variables | Description Value Reference 
M(0) Infants with maternal immunity at time t = 0 290 Assumed 
S(0) Number of susceptibles at time t — 0 800 Assumed 
vı (0) Vaccinated population taking one dose at time t — 0 300 Assumed 
V2(0) Vaccinated population taking second dose at time t = 0 | 150 Assumed 
V3(0) Vaccinated population taking third dose at time t = 0 100 Assumed 
I(0) Number of infected children at time t = 0 100 Assumed 
R(0) Number of recovered children at time t = 0 50 Assumed 
Parameters 
b Birth rate 0.0018 per day | [18] 
ш Natural death rate 0.0018 per day | [18] 
q Proportion of infants with maternal immunity 0.6 14] 
p Rate of departure from breast feeding 0.0111 2] 
B Contact transmission rate 0.00025 6] 
Un Reduction in risk from first vaccine 0.4 Assumed 
Na Reduction in risk from second vaccine 0.7 Assumed 
їз Reduction in risk from third vaccine 0.85 Assumed 
0 Recovery rate 0.1818 per day | [8, 15] 
Фф The daily rate of first dose vaccination 0.01667 Assumed 
T The daily rate of second dose of vaccination 0.00833 Assumed 
к The daily rate of third dose of vaccination 0.00556 Assumed 

Table 3b. Weight parameters of the Model 

Parameters | Description Value | Reference 

Bo Weight parameter for first dose 40 Assumed 

Bi Weight parameter for second dose | 40 Assumed 

Вә Weight parameter for third dose 40 Assumed 

Вз Weight parameter for treatment 5 Assumed 

















2.1. No control strategies. Under no strategy, Figures 2.1 and 2.2, 
shows the dynamics of rotavirus disease without treatment and vacci- 
nation. We consider a case where some of the children are breastfed 
while other are not. In Figure 2.1, breast-fed children move to class (a) 
and after a period of three months they join the susceptible class (b). 
Children that are not breastfed go directly to (b) class. At (b) children 
are vaccinated for first dose of rotavirus and they move to (c). We no- 
tice an initial increase in this class since they are coming from (b) and 
later a decrease because they are moving to class (d) for second dose of 
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vaccination. Again from (d) children are vaccinated for the third dose 
and they move to class (e) as seen in Figure 2.2. Again from Figure 2.2 
(f), we notice an increase in number of infected individuals who later 
decrease in number but after a long time, almost when t = 75 days. We 
lastly notice that children in the recovered class (g) take sometime to 
recover completely. This is due to the fact that, in absence of controls, 
the recovery level tends to be slow than in the presence of controls. 
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FIGURE 2.1. Profile (a) of the dynamics of rotavirus 
without controls 
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(e) (f) 
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FIGURE 2.2. Profile (b) of the dynamics of rotavirus 
without controls 


2.2. Optimal Treatment Only. Under this strategy, we optimize 
only treatment and all vaccination doses are put to zero, that is, ш = 
ио = из = 0. When treatment is optimized, we do not see any effect in 
the maternal class (Figure 2.3(a)-(1)), a slight increase in the number 
of susceptibles leaving this class but treatment has no direct effect on 
this class (Figure 2.3(a)-(2)). We notice an increase in the number 
of vaccinated children within Figure 2.3(a)-(3),(4), and Figure 2.3(b)- 
(5). Since treatment has a direct effect on the infected class (Figure 
2.3(b)-(6)) we notice that infection starts to decrease up to when t = 55 
days. This is well shown in Figure 2.3(c) where in absence of treatment 
infection takes longer to clear, that is, t = 75 days than in presence of 
treatment where it takes less time, t = 55 days. 
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Effect of treatment on the infected class 
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(c) Infected class with treatment only (red curve) and no treatment (blue curve) 


FIGURE 2.3. Plots (a) and (b) shows u4 profile where 
optimal control of treatment is indicated by the red curve 
and the blue shows no control while (c) shows the effect 
of treatment on the infected class (red curve) and no 
control (blue curve) 


2.3. Optimal Vaccination only. Under this strategy, in Figure 2.4, 
we optimize vaccination, considering all the three doses, that is, u1, ио, 
and из setting the treatment control variable u4 = 0. When vaccination 
is optimized, susceptibes take much longer in the susceptible class, that 
is, 50 days and in case of no vaccination, they take between 15 to 20 
days to move from this class. We further note that, vaccination has 
an indirect effect on the infected class, as shown in Figure 2.4(c), the 
disease disappear within the first 30 days, that is, at t = 30 and in 
absence of vaccination the disease takes 75 days to disappear. Again 
recovery takes 10 days when vaccination is implemented and without 
it, it takes 40 days for full recovery to take place. 
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(c) Infected class with vaccination only (green curve) and non- 
vaccinated (blue curve) 


FIGURE 2.4. A plot represents optimal vaccination only. 


2.4. Optimal Vaccination and Treatment. Under this strategy, in 
Figure 2.5 (a) and (b) we optimize both vaccination and treatment. 
When both treatment and vaccination are implemented, infection is 
controlled within the first 10 days as shown in Figure 2.5(c). 
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2.5. Optimal without control, with treatment only, vaccina- 
tion only, and both treatment and vaccination. In Figure 2.6, 
we deduce the effect of different control scenarios on the infected class 
since our aim is to reduce the level of infected class. Without controls, 
infection takes t=75 days for the infection to clear, with treatment 
only, t = 55 days, with vaccination only t = 30 — 40 days and in 
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Effect of vaccination and treatment on infected class 
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(c) Infected class with both treatment and vaccination (pink curve) 
and no vaccination, no treatment (blue curve) 


FIGURE 2.5. A plot represents optimal of both vaccina- 
tion and treatment. 


presence of both treatment and vaccination, t = 10 days. Thus we 
recommend both treatment and vaccination since infection disappears 
within a short time. 
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FIGURE 2.6. A plot represents optimal without control 
(blue curve), with treatment only(red curve), with vacci- 
nation only(green curve), and with both treatment and 
treatment (pink curve) 


3. CONCLUSION 


We derived and analyzed a mathematical model for the transmission 
dynamics of rotavirus disease that includes vaccination administered 
as three booster doses and treatment as our control strategies. We 
performed optimal control analysis of the model by applying the con- 
ditions necessary for optimal control of the disease derived using the 
Pontryagin's maximum principle. The numerical results show that in 
case of a disease outbreak we should used both control measures to 
reduce the infection burden within the first 10 days, otherwise, if we 
treat only, the infection will take 55 days to clear, if we vaccinate, the 
infection will take between 30 to 40 days to clear and if we do not 
control at all, the infection will take much longer time to clear, that 
is, 75 days. The later is very dangerous because if do not control the 
infection burden will increase hence making more children to become 
infected. 
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